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Q_r 1. Introduction 

Quantum Electrodynamics (QED) is one of the most successful and precise theory 
in physics. Already studied in the early years of quantum mechanics by Dirac, 
Pauli, Heisenberg, Fermi, Weisskopf and others, it was finally formulated in its 
definite form by Bethe, Dyson, Feynman, Schwinger and Tomonaga between 1946 
and 1950. This gave Feynman, Schwinger and Tomonaga a Nobel prize in 1965. 

The theory is the combination of quantum mechanics and Einstein's special 
relativity. It aims at describing the interactions between matter and light, at the 
microscopic scale where quantum effects are dominant. It allows to determine 
the time-dependent behavior of charged particles (like the electrons in an atom 
or a molecule) when they are coupled to photons (the quanta of light). Quantum 
Electrodynamics has an important symmetry called charge conjugation. The latter 
means that the behavior of positively charged and negatively charged particles is 
described in a similar manner. More importantly, the theory predicts that any 
charged particle automatically has an anti-particle with the same mass but an 
opposite charge (for the electron, this particle is called the positron), and that it 
is possible to create a particle/ anti-particle pair by providing a sufficient amount 
of energy to the vacuum. From Einstein's famous relation, this energy must be at 
least 2 x mc 2 . Because one can create matter from energy, the vacuum cannot be 
seen anymore as an empty and inert object as it is considered in everyday life. At 
the microscopic scale, the quantum vacuum is a fluctuating complicated system 
which participates to any physical phenomenon. 
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Quantum Electrodynamics is an extremely accurate theory. Its predictions are 
the most precise ever obtained from a physical model, when compared with experi- 
ments. The most famous successful predictions are the Lamb shift in the spectrum 
of the hydrogen atom and the electron anomalous magnetic dipole moment. The 
agreement with experiment is within a window of about 10~ 8 . 

In more technical terms, Quantum Electrodynamics is an abelian gauge field 
theory. The ideas of Dyson, Feynman, Schwinger and Tomonaga were later used 
to invent more complicated non-abelian gauge field (a.k.a. Yang-Mills) theories, 
like those describing the strong and weak forces. In units where the speed of light 
is c = 1 and Plank's constant is h — 1, QED only depends on two parameters 
whose value has to be determined by experiment: the mass m of the electron, and 
the Sommerfeld coupling constant a which is the square of the electron charge, 
a = e 2 . Of course, there might also be external fields which are applied to the 
system. The physical value of the coupling constant a is small (about 1/137) 
and, in the physical literature, Quantum Electrodynamics is always formulated 
as a perturbative theory. This means that the interesting physical quantities are 
formally expanded as a power series in a, and that only the coefficients of the series 
are computed explicitly, order by order. But there are divergences occurring and 
these coefficients are all infinite! In order to solve this problem, a regularization 
parameter A has to be introduced in the model. Then, the divergences are absorbed 
by making a change of variable for m and a, in order to obtain a well-defined formal 
series in the limit A — > oo. This procedure is called renormalization. 

It is fair to say that this perturbative formulation of QED is perfectly rigorous in 
the sense that there is a unique and well defined way to get the final answer (which, 
in most cases, is then in a surprising agreement with experiment). But, on the 
other hand, the perturbative nature is certainly frustrating from a mathematical 
point of view. In a famous quotation, Feynman himself said in 1985: 

The shell game that we play (...) is technically called 'renormaliza- 
tion'. But no matter how clever the word, it is still what I would call 
a dippy process! Having to resort to such hocus-pocus has prevented 
us from proving that the theory of quantum electrodynamics is math- 
ematically self-consistent. It's surprising that the theory still hasn't 
been proved self-consistent one way or the other by now; I suspect that 
renormalization is not mathematically legitimate. [H] 

There has been no dramatic change since Feynman's quotation, and renormaliza- 
tion has become a common (and somehow accepted) tool. 

On the mathematical side, there have been several works on models originating 
from QED, but not so many on the true theory itself. We should recall that one of 
the famous Millenium Prize Problems of the Clay Mathematics Institute concerns 
the construction of a well-defined Yang-Mills theory as well as the understanding 
of its low energy excitation spectrum. Yang-Mills theory gathers all non-abelian 
gauge field theories as alluded above, and any deeper understanding of the (abelian) 
theory of Quantum Electrodynamics is therefore desirable. 

In the last few years, mathematicians have been particularly interested in study- 
ing the interaction of light with non relativistic matter. In this simplified model 
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the electrons are quantum but they are assumed to have a speed which is much 
smaller than the speed of light, such that relativistic effects can be neglected. Their 
description then involves the Laplacian instead of the Dirac operator. Events like 
the creation of electron-positron pairs are encoded in the Dirac operator and they 
cannot be described within such a non-relativistic theory. But it is already a very 
important and fundamental problem to understand the effect of quantized light on 
non-relativistic particles. Some of the recent works concern for instance the exis- 
tence of a lower eigenvalue for the underlying Hamiltonian 7 ,25 ,34;, resonances 
and the relaxation to the ground state [212111] , problems related to divergences and 
mass renormalization (4,5 31,33 42 , and the stability of large systems [18.40. 43144) . 

Coming back to relativistic particles, some authors have considered Lattice 
QED [27115I]. Other studied linear and nonlinear models based on the Dirac op- 
erator for finitely many particles in an 'inert' vacuum and with a classical elec- 
tromagnetic field, see e.g. [TniHSlHIlUniHTJES] an d the references in [TB]. Finally, 
the quantum vacuum and the process of pair creation was investigated in a non- 
interacting setting (meaning without any light at all), e.g. in [36,37,45,461149]. 

With P. Gravejat, C. Hainzl, E. Sere and J. -P. Solovej, we followed another 
route in a series of works [^2 24,26 30; which originated from a fundamental paper 
of Chaix and Iracane [HUH] , and which was stimulated by the previous works [3"ll3"2"] . 
We considered relativistic particles in a fluctuating vacuum, both described by 
the Dirac equation, and in interaction with a classical electromagnetic field. That 
light is not quantized in this model prevents us from describing important physical 
effects. But, on the other hand, these seem to be the first mathematical results 
dealing with the quantum vacuum in interaction with light. Also, as we will 
explain, we are able to construct the associated model in a fully non-perturbative 
fashion, in a simplified setting. 

The purpose of this article is to review these recent results. We will particularly 
insist on the uncommon mathematical aspects, like those related to renormaliza- 
tion, and which are related to important physical effects. The reader interested in 
knowing more details can also read the last section of the review |16j . 



2. The quantum vacuum in classical electromagnetic fields 

In this section we present a model for quantum relativistic particles evolving in a 
fluctuating quantum vacuum, and interacting with a classical electromagnetic field. 
We will not discuss too much here the mathematical meaning of the equations, 
which will be explained in the next sections. Actually, most of the terms that we 
write here are infinite quantities if no special care is employed. We will concentrate 
on the time-independent model and look for stationary states. 

Our system is composed of particles (electrons or positrons) together with the 
quantum vacuum on the one side, and of a classical electromagnetic field {E 1 B) 
describing light on the other side, all evolving in the physical space R 3 . It is useful 
to use electromagnetic potentials V and A which are such that E = —W and 
B = V A A. We might also consider fixed external fields E cyLt = -V7 cx t and 



4 



Mathieu Lewin 



B C xt — V A A cxt which are applied to the system. The particles and the vacuum 
are described by the Dirac operator 

D m ,eA tot : =X^ a fc(~*7r eA to t(x)k^ + m/3 - eV to t{x), (1) 

fe=i Xk 

where A tot := (V to t,A tot ) with Vtot = V + V cxt and A tot = A + A cxt are the total 
electromagnetic fields, and e is the (bare) electron charge. The four Dirac matrices 
a. = (a>i, a>2, 0:3) and j3 are equal to 

ak '- = \a k "o) and ^ := (o -J 2 )' 
the Pauli matrices o~±, 01 and 03 being defined by 

^-(i 0)' CT2: =(° «') and CT3:= (o -1)' 

These matrices are chosen to ensure that (-D m ,o,o) 2 = — A + m 2 which is the 
quantum equivalent of Einstein's fundamental relation E 2 = c 2 p 2 + m 2 c 4 for the 
relativistic classical energy E in terms of the momentum p, with c = 1 in our case. 
The operator in (JTJ acts on the Hilbert space ft := L 2 (M 3 ,C 4 ). Under reasonable 
assumptions on Vtot an d .Atot, it is self-adjoint on the Sobolev space if 1 (R 3 ,C 4 ), 
see [53) . The charge conjugation is the anti- unitary operator defined on fj by 
c tof := i(3ct2f, and which satisfies ^Dm^A^^^ 1 = —D, n _ e A tat - This relation 
implies that the spectrum of -D m ,o,o is symmetric with respect to 0, 

ff(flm,0,o) = (-00, -TO] U [TO, 00). 

When Vtot and A to t decay at infinity, the essential spectrum stays the same by the 
Rellich-Kato theorem, o~ ess (D mie A tot ) — (—00,— to] U [to, 00). Isolated eigenvalues 
of finite multiplicity can however appear in the gap (—to, m) when Vtot 7^ 0. 

The state of one particle is described by a normalized wave function ip € 
L 2 (M 3 ,C 4 ), f R3 \ip\ 2 = 1, and its corresponding energy is (<p, D m ^ e A tot <f) (which 
has to be understood in the form sense when <p e H^ 2 {M. 3 ,C 4 )). The fact that the 
Dirac operator is unbounded from below, hence that the energy can be arbitrarily 
negative was very surprising for his inventor, P. A.M. Dirac [101112) . In quantum 
mechanics it is always assumed that the most stable state of the system is the one 
with the lowest energy, and there does not seem to be any here. From a mathe- 
matical point of view this makes a huge difference as compared to non-rclativistic 
models based on the Laplace operator —A. That the energy is unbounded from 
below means that minimization methods cannot be employed and that one has to 
resort to complicated min-max techniques to construct solutions [16] . 

Dirac did not want to renounce the physical picture that states with lower 
energy are more stable. So, in 1930 he suggested to reinterpret the problem by 
changing the role of the vacuum: 
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We make the assumption that, in the world as we know it, nearly 
all the states of negative energy for the electrons are occupied, with 
just one electron in each state, and that a uniform filling of all the 
negative-energy states is completely unobservable to us. [IT] 

Physically, one therefore has to imagine that the vacuum (called the Dirac sea) is 
filled with infinitely many virtual particles occupying the negative energy states. 
With this conjecture, a real electron cannot be in a negative state. This is because 
electrons are fermions and that two such particles can never be in the same quan- 
tum state (Pauli principle). With this interpretation, Dirac was able to conjecture 
the existence of "holes" in the vacuum, interpreted as anti- electrons or positrons, 
and which were later experimentally discovered by Anderson [2]. Dirac also pre- 
dicted the phenomenon of vacuum polarization: In an external electric field, the 
virtual electrons are displaced, and the vacuum acquires a non constant density of 
charge. The idea of a fluctuating quantum vacuum was born. 

Let us now explain Dirac's idea in more mathematical terms. The state of N 
electrons can be represented by N wave functions (pi, —,<Pn G i? 1 (R 3 ,C 4 ), such 
that (ipi,ipj) L 2 — 5ij, the latter constraint being the mathematical formulation of 
the Pauli principle. The total energy is the sum of the energies of the individual 
electrons, which can be written as 

N 

{Vi,D m ,eA M ^3) = Tr (D m ,eA toi P), 

where P is the orthogonal projector onto the A^-dimensional space spanned by the 
V?j's. Using the bra-ket notation we have P = 53j=i IViXVjl- Now if we allow any 
number N of particles and look for the state of lowest energy, the formal solution 
is the negative spectral projector 

P = P rn,eA tot ■= l(-oo,0) (^Mtot) • (2) 

Think of a finite hermitian matrix M, then the solution to the minimization 
problem inf{Tr(MP), P 2 = P = P*} is the negative spectral projector P = 
l(-oo,o) °f M, and the corresponding energy is — Tr M_ with X- := — min(x, 0) 
denoting the negative part. The minimizer is unique when ker(M) = {0}. In our 
case the state in @ corresponds to filling with particles all the negative energies 
of P> m ,eA tot as suggested by Dirac. This state P~ eAtat 1S the quantum vacuum in 
the presence of the electromagnetic potential A to t ■ It has infinitely many particles 
and the corresponding total energy is also infinite. 

Remark 2.1. There are several ways to give a mathematical meaning to the 
assertion that P~ eA t minimizes P H> Tx{D m , e A tot P). The first is to use a so- 
called thermodynamic limit. The ambient space fj = L 2 (R 3 ,C 4 ) is approximated 
by a sequence of finite-dimensional spaces in which everything makes sense and 
it is proved that in the limit the minimizer converges to P~ eA . The second 
possibility is to directly argue in the whole space fj that 

TrP m , eAtot (P-P meAtot )>0 
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for all projection P. Under appropriate assumptions on P (for instance P — P m 
finite rank and smooth), we can write 

Tr An : eA tot (P — P m ,eA tat ) = Tr I An,eA tot | {Pm,eA tot (P ~ ^m,eA to ,)^m,eyl tot 

_ P- (P — P~ )P~ ] 

m,eA tot V m.eAtot'' m,eA tot J 

= Tr |-D m , e A tot r^^-Pm.eAtot-f-fm.eAtot 

+ ^,eA tot (1 ~ ^,eA tot )l A^eAj 172 > 0. (3) 



We have used here the commutativity of the trace and that D mje A tot — \Dm,eA tot \ 
(-^m e A e A ) ' by definition of the spectral projections. The last term in ([3]) 

is the trace of a non- negative operator which is always well-defined in [0, oo], and 
which can be taken as a definition for the trace on the left side. See [2l)l Sec. 2.1] 
for a systematic and abstract theory of traces of this form. 

When A to t = 0, then P~ is nothing but the negative spectral projector of 
the free Dirac operator P m ,o, which represents the free vacuum. This operator is 
translation invariant and, even if the corresponding charge density is infinite, it is 
somehow constant. This is the 'uniformity' alluded to in Dirac's quotation. When 
eAtot 7^ 0, the state of the vacuum changes. Later we will want to optimize over 
A with A ext fixed, but before we have to modify a bit our theory in order to make 
it charge-conjugation invariant. If we take an arbitrary state P = P 2 , then we get 
Tr {D m e A tat P) = — Tr (^D m ^ e A tat '^P'^'~ 1 )- This is not very satisfactory because 
— ^P^ -1 is not a fermionic state and we cannot reinterpret this as the energy 
of something. It is better to subtract half the identity to the operator P, which 
amounts to adding a (infinite) constant. So we consider instead the energy 

TrD m ^ Atot I P- - I = TrD m ,eA tot ^ 

withP- 1 := I-P. Now we see that Tr D m ^ Atot (P- P^) = Tr D m ,- e A tot {P'-P'^) 
with P' := ^P^- 1 . The subtraction of half the identity to P is a common 
technique for systems at half filling like ours. It is explained using the formalism 
of second quantization in [50] . 

We are now able to write the total Lagrangian of our system, which is the 
(formal) sum of the particle energy and Maxwell's classical Lagrangian: 

jSf m , e (P,A;A ext ) :=Tr^„ i!e(A+Aoxt) (p - + -L \B\ 2 - |P| 2 , 

where we recall that E = — VV and B = V A A. Our purpose is to look for critical 
points of this Lagrangian, obtained by minimizing over P and A and maximizing 
over V, as is usually done in classical electrodynamics. We have already explained 
that the minimum over P with A fixed gives P = P^ eA . It is also possible to 
optimize over A, with P fixed. The optimal potentials solve Gauss' equation 

-Ay = 47repp_i/2, -AA = An e j P -i/ 2 
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where Pp-1/2 an d jp-1/2 are the (formal) density of charge and density of current 
of the vacuum in the state P. These are defined by Pm{x) — Tr C 4 M(x, x) and 
jM(x)k — Tr£4 a.kM(x, x), for any locally trace-class operator M with integral 
kernel M{x, y). Note that these are zero for the free vacuum: 

as is shown using the commutation properties of the matrices a^. This emphasizes 
the importance of subtracting half the identity to our vacuum state. 

Our task is to optimize over both P and A and it is reasonable to think that 
the free vacuum P = P~ is the optimal state when A cxt — 0. But, the functional 
Jzf m . e is infinite and so we cannot easily optimize it. We can however subtract the 
universal (infinite) constant «5f mje (P~ , ; 0) and consider the relative Lagrangian 

#£ e (P, A ; A cxt ) = Tr D m>0 (P - P mfl ) - e [ j p _ p - ■ (A + A cxt ) 

-el P P -p-(V + V cxt ) + ±[ \Bf-\Ef. (4) 

JR3 m '° ™ JR3 

Because the difference P — P~ is a much better behaved operator than P — 1/2, 
we will be able to give a clear meaning to this Lagrangian. The purely electrostatic 
case A = A ext = is much easier to deal with than the general case, and we discuss 
it first in the next section. 



3. The purely electrostatic case 

In the relative Lagrangian we now assume A = A cx t = and we maximize 
with respect to V. We find the following energy functional 

f a 
TrD mi oQ-e/ PQV cyLt + -D{p Q , p Q ) 

JR3 1 

with a = e 2 (the Sommerfeld coupling constant), Q := P — P^ and 



Jr3 Jr3 \x - y\ J R 3 \k\ z 

which is the Coulomb (i.e. H~ 1 (M 3 )) scalar product. The associated space will be 
denoted by C :— {/ G S*' (K 3 ) : D(f,f)< 00}. In our setting we are interested in 
taking for V ey ± the external field induced by some localized density of charge v ext , 
describing for instance the nuclei of a molecule. So we take V cx t = ev cx t * \x\ _1 
and get the energy 

C°,a(Q) = Tr D m ,oQ - aD(p Q ,u ext ) + ^D(p Q , PQ ) (5) 



which is called the reduced Bogoliubov-Dirac-Fock energy [81127] . 
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Under suitable assumptions on v cx t, the minimization of this functional in Q 
makes perfect sense. First we remark that, by the Cauchy-Schwarz inequality, 

OL (X 

-aD(p Q ,v ext ) + -^D(p Q ,p Q ) > --D(v cxt ,v cxt ) 

which is finite when v cxt € C (for instance for v cxt & L 6 / 5 (R 3 ), by the Hardy- 
Littlewood-Sobolev inequality [41]). Then we recall that TrD m oQ > for any 
operator Q = P — P^ , provided the trace is defined like in ([3]). This is reminiscent 
of the fact that P~ formally minimizes P n- Tr D m ^P. From all this we conclude 
that £ UBXt (Q) > — aD(v cxt , ^ext)/2 and the energy is bounded from below. 

However we are unlucky that the infimum of £^ is never attained [271 Thm 
2], except when v cxt = where Q = is the optimizer. The clear mathematical 
statement is 

inf (TrD mfi Q-aD(p Q ,u) + ^D(p Q ,p Q )) = ~D{v, v). 

-P~. <Q< P ,t,0 I i 

In the infimum one can restrict to finite-rank smooth operators Q, which ensures 
that all the terms make sense. We see that there is never a minimizer when v cxt ^ 0, 
because it should satisfy both pq = v and Q — 0. As explained in [27] , this issue 
is due to certain divergences in Fourier space at infinity. 

In order to give a meaning to our minimization problem, we have to impose a 
cut-off for the large Fourier frequencies. There are several ways to do this, some 
being more natural than others. In the purely electric case, we can consider a 
simple sharp cut-off consisting in replacing our ambient Hilbert space f) by 

ru:={/ei 2 (R 3 ,C 4 ) : f(k) = for |fc| > a} . 

Note that D m ^ maps $j\ into itself. For simplicity we use the same notation D m ^ 
and P^ for the associated restrictions to The operator Q is defined on Sj\ 
and its energy has the same expression as in ([5]). It now has minimizers. 

Theorem 3.1 (Existence of the polarized vacuum [2l)ll2"7] ) . Let m > 0, 0<A<oo 
and a > 0. For any fixed i/ cxt G C = H (M. ), the minimization problem 

inf {£%%(Q) : Q : io A -> «a, -P m ,o <Q< i£,o> T Tr < oQ<o < ™} (6) 

admits at least one minimizer . All these minimizers share the same density 
PQ r . Any minimizer is a solution of the nonlinear equation 

= l(_oo,o) (D m )-P- >0 +S 
D* = Ii A (^D mfi + a(p Qr - Vext) * |x| _1 )Ha 

where II a is the projection onto $)a, and < 5 < 1{ j (-D*). 

// an 1 / 6 2 11 / 6 D{is,v) 1 / 2 < y/m, thenkei(D^) = {0} and Tt(P+ Q*P+ + 
P~ Q*P m ) = 0. Hence S = in ([Jj, and the minimizer Q* is unique. 



(7) 
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The equation ([7]) on the infinite rank operator P* = + P~ can be inter- 
preted as an infinite system of coupled nonlinear partial differential equations. As 
we have said we always use the convention 

TrD m ,o<3 := Tr |£> m ,o|(Q ++ - Q ) 

where Q++ := P+ fi QP+. > and Q— := P m , QP m , < when Q satisfies 
the constraint -? m0 < Q < P^ l0 - Since with the cut-off A the operator D m o 
is bounded, it is sufficient to ask that Q ++ and Q are trace-class in order to 
properly define the trace, as required in ([5]). On the other hand simple algebraic 
manipulations show that the constraint is equivalent to 

- P-, <Q< P+ o Q 2 < Q ++ - Q— (8) 

as remarked in [3]. So we see that we must also have TrQ 2 < oo. In other 
words, for the first term in the energy to be finite, Q must be a Hilbert-Schmidt 
operator whose diagonal blocks Q ±±: are both trace-class. Under these conditions 
and with the ultraviolet cut-off A, it was proved in [23 that pq 6 L 2 (K 3 )nC. This 
is sufficient to prove Theorem 13.11 using simple arguments from convex analysis. 
In [26l[27] a more complicated model with an additional so-called exchange term 
is considered and this term makes the analysis much more involved. 

Let us emphasize that Theorem l3.1l is valid for any value of the coupling constant 
a and any value of the ultraviolet cut-off A. This is therefore a non-perturbative 
result. Because no smallness assumption is made on the external density z/ ex tj 
the model is appropriate for the description of non-perturbative events like the 
spontaneous creation of electron-positron pairs in strong external fields [50] . Note 
that for simplicity we have minimized the energy with the relaxed constraint ([8"jl 
instead of the original constraint that P := Q + P m = P 2 . However, for a small 
external density f ex t> the minimizer is unique and the corresponding state P* 
is an orthogonal projection, as expected. 

In Theorem 13.11 we have constructed the polarized quantum vacuum in the 
presence of the external density f e xti which is by definition a global minimizer of 
the energy. It is also possible to construct states having a certain number N of 
'real' electrons (or positrons). More precisely, we can minimize the energy under 
a charge constraint of the form "TrQ" := Tr(Q++ + Q — ) = N, see [231128] , 
Minimizers do not always exist, depending on the value of TV and on the strength 
of the density v e xt- Estimates on the maximal number of electrons that can be 
bound by ^ ox t were provided in |23| , following ideas of Lieb |39] . Any minimizer 
with charge constraint solves a nonlinear equation similar to Q, with the energy 
level replaced by a Lagrange multiplier /i G {—m, to): 

Q, = l ( _ OOi/i) (D.)-P- + * (9) 

where now < 6 < 1 {ai} (£>*). If /Lt > 0, splitting l(-oo >A i) = 1 (-oo,o) + l(o, M ) S ives 
the corresponding states for the polarized vacuum and for the electrons. If /i < 0, 
we split l(_oo,ju) = l(-oo.o) ~ 1(m,o) an d get the state of the polarized vacuum 
and of the positrons. In the rest of this article we will not discuss further this 



10 



Mathieu Lewin 



constrained minimization problem, and we will instead concentrate on the pure 
vacuum case considered in Theorem 13. II 

The minimizers obtained in Theorem 13.11 are very singular mathematical ob- 
jects. The following says that is never a trace-class operator. 

Theorem 3.2 (Renormalized charge |23j). Let Q* be a minimizer as in Theo- 
rem lff.il for some v cxt € L^R 3 ) DC. Then we have pq, G L X (R 3 ) and 



MQ++ + Q--)- [ v cxt 
Jm. 3 



PQ, ~ ^cxt = ——= ^ (10) 

1 + atSx/m 



where 



1 f ^ m 2 +A 2 z 2 -z i /i J 2 A 5 21og2 „..«. , . 

B\/ m = - / -r: 5-^ = ^- log - - — + -^2- + O m 2 /A 2 . 11 

7r 7 1 — z z Sir m 9tt Sit 

In particular, if a r K 1 ^2 11 ^D{v,v) 1 ^ 2 < \fm, then is not trace-class. 

Let us recall that a compact operator B is trace-class when Tr \B\ — Tr %/ B*B < 
oo. This is equivalent to saying that J2j (Vji l-^IVi) < 00 f° r one (hence for all) 
orthonormal basis {ipj} of the ambient Hilbert space. By the spectral theorem, 
this is also the same as saying that the eigenvalues of B are summable, ^ \Xj\ < 
oo. If the ambient space is an L 2 space, we can write by the spectral theorem 
B = J2j ^j\'4 ! j)(i } j\ an d we get that the corresponding density ps = X)j Aj'IV'jl 2 is 
in L 1 and that J pb = TrB. 

For a non positive operator B, it can happen that Y] ^ (ifj, Bipj) is convergent 
for one basis and not for another. This is exactly what is happening for our 
operator Q*. The two diagonal blocks Qf + and Q~~ are trace-class, which means 
that J2j { < Pj' , ^*f > t) (^7' Q*^Pj) converges for any orthonormal basis {<pf} of 
-f)^ 1 := P^ $j\. The surprise is that pq t is always in L 1 (R 3 ) but that in general 
Jr3 PQ, Tr(Q ++ + Q ). The discrepancy between these two quantities is 
universal and it is given by the relation (|10p . The problem comes from the off- 
diagonal densities Pq±t which are in L 1 (R 3 ) but do not have a vanishing integral. 
More precisely, only the first order term when Pq±=f are expanded in a power series 
in a contribute to J R3 pq t . The proof of Theorem consists in studying this term 
in details [23] . 

Theorem 13.21 has a natural interpretation in terms of charge renormalization. 
Imagine that we put a nucleus of charge Z = L 3 f cxt in the vacuum, which is weak 
in the sense that D(v ext , f cx t) is small, and let P* = + P~ be the correspond- 
ing unique polarized vacuum obtained by Theorem 13.11 This vacuum is neutral, 
Tr(Q^ + + ) = 0, which means that the external field is not strong enough 
to create electron-positron pairs. In reality we never measure the charge of the 
nucleus alone, but we always also observe the corresponding vacuum polarization. 
Hence we do not see Z, but rather Z/(l + aB^/ m ). This corresponds to having a 
physical coupling constant a p h given by the renormalization formula 

a + a P h tir,\ 

«ph = -, , „ ^a=- . (12) 

1 + a±S\/ m 1 - Ckph-DA/m 
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In our theory we must take a p h — 1/137, the 'bare' a can never been observed. 
So using the change of variable (fTiZj) we should express any other physical quantity 
predicted by our model in terms of ev p h, m and A only. 

The natural question arises whether it is possible to remove the ultraviolet cut- 
off A, by keeping a p h and m fixed (note that the mass m is not renormalized in this 
model). The answer is clearly no! We immediately get from (TT2j) that a p h,B^/ m < 1 
and therefore a p h — > whenever we try to take A — > oo. This phenomenon is called 
the Landau pole |38j and one has to look for a weaker definition of renormalizability. 
The cut-off A which was first introduced as a mathematical trick to regularize the 
model has actually a physical meaning. A natural scale occurs beyond which the 
model does not make sense. Fortunately, this corresponds to momenta of the order 
A ~ me 37r / 2Qph , a huge number for a p h — 1/137. 

In [23], the regime a p h < 1, A > 1 with a p hlogA fixed was studied. It was 
proved that the dressed density of the nucleus admits a Taylor expansion 

K 

a(v ext ~ PQ,) = "ph^cxt + ^2 v k a^ 1 + 0(a^ 2 ) (13) 

k=l 

where the terms v\. in the expansion are independent of the value of a p h log A. 
The first correction v\ gives rise to the famous Uehling potential [SUES]. The 
relation (|13p shows that the density of the nucleus can be renormalized order by 
order. It is believed that the series in (TT3)) is divergent [14], but no mathematical 
argument has been provided so far. 

To summarize, in the purely electrostatic case we have a well defined non- 
perturbative theory for all m > 0, all a p h > and all cut-off A such that 
a p hi?A/m < 1- In the regime a p h « 1, A » 1 with a p hlogA fixed, the per- 
turbation series computed by the physicists is recovered. 

4. The electromagnetic case 

Let us go back to our relative Lagrangian ((4]), including now the electromagnetic 
fields i?tot = B + B CKt - Optimizing with respect to P we find a formal Lagrangian 
depending only on the classical fields 

J2f£UA ; A«t) = \Tr (\D mfi \ - \D mMA+A ^ t) |) + -L J \B\ 2 - \E\ 2 . (14) 

This Lagrangian is not really well-defined because of the ultraviolet divergences 
which have to be taken care of. On the contrary to the pure electrostatic case 
studied in the last section, we cannot impose a sharp cut-off here, as it is extremely 
important to keep the magnetic gauge invariance corresponding to replacing A by 
A+Vip. In [22] we used the Pauli- Villars method [15] > which consists in introducing 
two fictitious particle fields of very high masses mi, 3> 1 which play the role of 
ultraviolet cut-offs. The Pauli- Villars-regulated Lagrangian is defined by 

2 

if^ (A;Aext) = 1 Tr ^ c .^| Dm . |_| jDm . ie(A+>le5tt) |) + i.^ \B\ 2 -\E\ 2 (15) 
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where cq = 1, mo — m and the fictitious fields are chosen such that 

2 2 

E c i=E c X = °- ( 16 ) 

The main result of [25] is the following 

Theorem 4.1 (The Pauli-Villars- regulated vacuum in electromagnetic fields 22 ). 
Let e 6 1, m > and Ci,C2,TOi,m2 satisfying (I16[) . 

(i) XTie functional A H- Tr Trc4 X)j=o c i (J-^TOj.oI — |-Dmj,A|J * s well-defined for 
A G C^°(R 3 ,R 4 ) and it admits a unique continuous extension to 

H^ iv (R 3 ) — {(V,A) e L 6 (M 3 ,M 4 ) : divA = 0, (-VV, V A A) € L 2 (R 3 ,R 6 )} . 

So the Pauli- Villars-regulated Lagrangian (|15p is well-defined and continuous on 



(ii) For any A cxt £ ij^ iv (R 3 ) with e\\A ext \\ ^i < r^frn/2, there exists a unique 
solution A* = (14, A*) € iJ^ iv (R 3 ) to the min-max problem 

J?™(A*,A cxt )= max inf ^ v e (A;A ext ) 

||VV||r,2< r ^ \\VAA\\ L 2<?^ 

= min sup ££^{A ; A ex t)- 

||VAA|| X , 2 <^ ||vV|| i2 <^ 

/or some radius r > which only depends on Ylj=o l c il( TO / m i) an< ^ provided that 
e stays in a bounded set. 

(ii) The four-potential A* = (V*, A*) is a solution to the nonlinear equations 

{ -^ = 4 / e ^.+^ t , (18) 

\ -AA. = 47rejA.+A cxt , V ' 

where PA„+A cxt and JA„+A cxt re/er to t/ie charge and current densities of 

2 

Q* = T^Cj l(-oo,0)(An J -,e(A»+A oxt ))- (19) 
3=0 

(iii) If A cxt = 0, then A* = as well. 

This result completely settles the problem to construct Dirac's vacuum in small 
external fields. It is an open question to do the same in large external fields (which 
could be useful to understand non-perturbative events like the creation of electron- 
positron pairs). It is also not known whether our solution is global. For the free 
vacuum, A cxt = 0, we conjecture that the saddle point A = is the unique 
solution to ([TT]) with r = oo. Even if this was not done in [22J, it is possible to 
renormalize the coupling constant a in the same fashion as in the previous section. 
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